High Fidelity Adiabatic Quantum Computation via Dynamical Decoupling 
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We introduce high-order dynamical decoupling strategies for open system adiabatic quantum computation. 
Our numerical results demonstrate that a judicious choice of high-order dynamical decoupling method, in con- 
junction with an encoding which allows computation to proceed alongside decoupling, can dramatically enhance 
the fidelity of adiabatic quantum computation in spite of decoherence. 



Introduction. — In adiabatic quantum computation (AQC) a 
problem is solved by evolving in the ground state manifold of 
an adiabatic Hamiltonian H^(t), with t e [0,T] |H1|2). Th e 
ground state of the beginning Hamiltonian Hp = -ff a d(0) is 
assumed to be easy to prepare, while the final one, the ground 
state of the problem Hamiltonian Hp = iJ at j(T), represents 
the solution to the computational problem. AQC has been 
shown to be computationally equivalent to the standard cir- 
cuit model of QC OJd, and is being pursued experimentally 
using superconducting flux qubits [8] and nuclear magnetic 
resonance J9JQID- However, in spite of evidence of intrin- 
sic robustness ifTTlTTBI and proposals to protect AQC against 
decoherence ifTTl [T8l . AQC still lacks a complete theory of 
fault-tolerance, unlike the circuit model of QC [ 1 9—24 1 . In 
fact, even identifying an acceptable notion of fault-tolerant 
AQC (FTAQC) is an open problem. In the circuit model a 
fault-tolerant simulation allows one to generate the output of a 
given ideal circuit, to arbitrary accuracy, using the faulty com- 
ponents of another circuit [25 1. A similar definition for AQC 
would presumably involve the simulation of an ideal adiabatic 
evolution using a faulty one, but if "faulty" is simply taken to 
mean "non-adiabatic", then one can just use the equivalence 
proof [3-7 1 and circuit-model fault-tolerance [19-24] to ar- 
gue that the problem is already solved. However, this argu- 
ment misses the point since to qualify as AQC, at least the 
computation should remain adiabatic, i.e., the defining fea- 
ture of AQC — the adiabatic preparation of the ground state of 
Hp — should be preserved. On the other hand, it seems too re- 
strictive to require the techniques used to address decoherence 
and noise to be adiabatic as well. We thus propose the follow- 
ing characterization of FTAQC: 'Given a closed-system AQC 
specified by Hb and Hp, and e > 0, a fault tolerant open- 
system simulation will use adiabatic evolution between two 
faulty, encoded Hamiltonians Hp and Hp derived from Hp 
and Hp, so that the final system-only state of the simulation is 
efficiently decodable to a state that is e-close (in fidelity) to the 
ground state of Hp. In addition, the simulation may involve 
any other faulty non-adiabatic error-correction, suppression, 
or avoidance operations.' 

Our characterization is meant to convey that the computa- 



tion should be adiabatic, and apart from that the error cor- 
rection can be anything. We make no attempt to rigorously 
quantify the relation between the "ideal" and "faulty" pairs 
Hp, Hp and Hp, Hp, the types of allowed error-correction 
operations, or the nature of the decoding step. Instead, we 
demonstrate here that it is possible to approach ideal AQC in 
an open system, using ideas guided by our characterization 
of FTAQC. Specifically, we shall assume that the underlying 
computation is indeed adiabatic but encoded into new Hamil- 
tonians, and that the protection is non-adiabatic. 

Ref. [ 18 1 introduced an AQC protection method which fits 
this FTAQC approach. Protection is carried out by means of 
dynamical decoupling (DD) [26-28], a non-adiabatic tech- 
nique which involves the application of strong and frequent 
pulses to the system, designed to decouple it from the envi- 
ronment, or bath. To ensure the compatibility of DD with 
AQC, all qubits are encoded into a quantum error detecting 
code 11291 , which allows DD pulses to be applied that com- 
mute with the adiabatic evolution, while at the same time act- 
ing to decouple the system from the bath [30]. Ref. [ 18] relied 
on first order DD sequences, resulting in a tradeoff between fi- 
delity and DD sequence bandwidth, and conjectured that high- 
order DD, in particular concatenated DD (CDD) [31 1, should 
alleviate this tradeoff. 

Here we demonstrate, using numerical simulations, that fol- 
lowing the strategy of Ref. [ 18 1, but using high-order DD se- 
quences, and in particular CDD, it is possible to dramatically 
enhance the fidelity of AQC in an open system setting. Our re- 
sults support the idea that FTAQC, in the sense characterized 
above, is indeed attainable using appropriately chosen DD se- 
quences. 

Algorithms. — We consider two well-known AQC algo- 
rithms, which correspond to first and second order quantum 
phase transitions, respectively, and thus to different challenges 
for AQC IB] [321 133). The first is Graver's algorithm for the 
identification of a marked element in an unsorted list of N 
elements, using the minimum number of oracle queries l34l . 
This can be done in 0(s/N) queries, which is a quadratic 
improvement over the best possible classical algorithm J35). 
Recast in the language of AQC l36l[37l . Graver's algorithm is 
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denned by the n-qubit Hamiltonian H^(t) = [1 - f(t)](I - 
\u)(u\) + f(t)(I - \m)(rn\), where \u) denotes the uniform 
superposition over all N = 2™ computational basis states, \m) 
is the marked state, and / is the identity operator. The mini- 
mum spectral gap A^ in = 0(l/\/~N), and the total run time 
required to reach the ground state |$q (T)) = |m) of H^ d (T) 
is T ~ 0(\ZN) provided the optimized interpolation function 
f(t) = \ - 2 ^L_ tan[(l - 2t/T)axccos(l/vOV)] is used 

ESI El. 

The second algorithm we consider is 2-SAT on a ring [1|. 
Given a set of clauses {Cj}" =1 associated with n bits, such 
that each clause acts only on adjacent bits, the 2-SAT Hamil- 
tonian acquires minimum energy when there are satisfying 
assignments for all clauses. In the case that the clauses de- 
fine agreement between adjacent bits (00 or 11, but not 01 
or 10), the 2-SAT Hamiltonian is that of the transverse-field 
Ising model with periodic boundary conditions af> +1 = of: 

H*l AT (t) = (l - |) - of) + £ £™=i W - 

Here oj and erj are the spin- 1/2 Pauli matrices acting on qubit 
j. The symmetric ground state of H^ AT (T), |$2 SAT (T)) = 
(|0 . . . 0) + |1 . . . l})/%/2, satisfies the clause agreement con- 
dition. The minimum gap occurs near t = 2T/3 and A m ; n = 
0(l/n); it can be found using the standard fermionization 
method flT] (38j. The linear interpolation in H^ AJ (t) is not 
optimal; see Ref. [39| for the optimal path. The 2-SAT on a 
ring problem is not associated with a quantum speedup (there 
exists a classical linear-time algorithm for 2-SAT (40)), but 
is instructive nonetheless, and important in its own right as it 
corresponds to the preparation of the maximally entangled cat 
state \$f AT (T)}. 

Error model. — To keep our numerical simulations manage- 
able, we considered interactions with a classical environment. 
Thus the "faulty Hamiltonian" is 

H Q (t)=H ad (t) + H en (t), (1) 

where the error Hamiltonian comprises interactions between 
the qubits and time-dependent classical fields: 

n 

E E^K- ( 2 ) 

fte{x,y,z} j=l 

Each field £j(t) is a zero-mean random Gaussian process with 
spectral density S^{ui) = 5jk5^ v S{u) [41 ], where 

S(u>) = -i= [ °° (*)e£ (* + r))e*"dr = ~j=e^ )\ 

(3) 

where (•) denotes a Gaussian ensemble average. Note that (3 
plays the role of the spectral cutoff (1//3 is the bath correlation 
time). While H sn (t) is not as general as a quantum environ- 
ment, it still represents an interesting and relevant error model, 
e.g., due to charge noise in superconducting qubits l42l l43l : 
its terms combine with H a <i(t) to produce faulty operations by 
inducing random time-dependent fluctuations in the spectrum 
of H a i(t). In particular, these fluctuations also act to modify 
the size and location of the minimum energy gap A m [ n of _ff a d. 



Dynamical decoupling. — The problem we attempt to solve 
using DD is to perform high fidelity AQC in spite of the pres- 
ence of H en . The decoupling pulses are introduced through an 
additional time-dependent control Hamiltonian Hc(t), which 
generates a unitary pulse propagator Uc(t). We consider 
zero-width pulses separated by finite intervals. The inclu- 
sion of pulse-width errors is left for a future study focusing 
on a more complete picture of fault tolerance; CDD is known 
to be relatively robust against such errors QT1 l44l [451 . The 
Hamiltonian H(t) = Hc(t) + H (t) generates the complete 
dynamics of the system in the presence of DD, represented by 
the unitary evolution operator U(t,0). To suppress H en: (t), 
while preserving H^t), we require that each term in H en -(t) 
anticommute with some pulse operator comprising Hc(t), 
while [Hc(t), H a &(t')] = Vt, t'. Upon satisfying these con- 
ditions the time evolution operator in the interaction ("tog- 
gling") frame with respect to Hc(t) is 

U(T,0):=U f c (T)U(T,0) (4) 

= T e -iT Jo H*{s)d S + 0[(\\H e jT) a+1 ], 

where T denotes time ordering, and H elT is an effective error 
Hamiltonian, which can be computed using the Magnus or 
Dyson series Poll . DD becomes effective provided the "noise 
strength" ||ij err ||T < 1 fl4ll . The larger the "decoupling order" 
a, the closer to ideal is the adiabatic evolution. Previously 
only the case a = 1 was analyzed ifTSIl ; we shall show that 
protection via high-order DD sequences, in particular CDD, 
alleviates this tradeoff. 

Stabilizer decoupling. — As in Ref. [ 18 1, to satisfy the non- 
interference condition [Hc(t), H Rl i(t')] = we make use of 
the [[n,n - 2,2]] stabilizer code C, encoding n - 2 logical 
qubits (n even) into n physical qubits [29 , 47 1. The stabilizer 
of C is S = {I,X,Y,Z}, where X(Y,Z) = ®" =1 af v,z) . 
The encoded single-qubit operators are ctJ = cr^crj +1 and 
ct| = <Tj +l a*, where j = 1,2, ...,n- 2. AQC over C is 
implemented by replacing each Pauli matrix in H. d i{t) by its 
encoded version, yielding an encoded adiabatic Hamiltonian 
H a d{t) which is fully 2-local. 

CDD for AQC. — Consider a unitary decoupling group Q = 
{#fc}fcli> chosen to implement first-order decoupling (a = 1). 
DD effectively averages out Herr by symmetrizing it over 
Q Gil EH). We choose Q = S and hence G = 4 for the 
[[n,n - 2,2]] code. CDD achieves higher order decou- 
pling by concatenating this symmetrization, so that each ad- 
ditional level averages out the remaining leading order. So far 
CDD has been defined only for piecewise constant Hamilto- 
nians EU |44l |48]|. Adapting CDD to incorporate the time- 
dependence of H(t), we have at the I + 1th level of concate- 
nation 

Q 

4dd( t .o)= Us^gl, i>0- (5) 

fcl=l 
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FIG. 1. (color online) Trace-norm distance between the CDD- 
protected final state p(T) and the desired ground state \<&o(T)), as 
a function of total run time T for Grover's algorithm, in units of the 
inverse minimum gap. Cutoff frequency: /3 = A/5. The ideal (solid 
black) and faulty (empty squares) evolutions are included for refer- 
ence. Insert shows a close-up for large T, using a log scale for the 
vertical axis. Performance improves monotonically with concatena- 
tion level I, with the corresponding sequence denoted CDD; . Error 
bars are due to averaging over 30 random realizations of H m (t). 
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FIG. 2. (color online) As in Fig. [TJ for the 2-SAT on a ring prob- 
lem. The curves for ideal evolution and CDD4 overlap to within 
our numerical accuracy up to T = 25/ A^^ T . Note the minimum 
at T « 12A„fj AT for the faulty evolution, suggesting the existence 
of an optimal open system evolution time. This optimal time in- 
creases with concatenation level, until it disappears in the ideal case 
and for CDD4. CDD boosts the deviation by a factor of 10 from 
D « 0.2 (at T ~ 12.5/A^f n AT ) in the faulty case to D « 0.02 (at 
T ~ 34.5/A^ n AT ) for CDD 4 . 



Starting from m = 2, is calculated recursively from 
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with DD-free evolution segments 
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where the DD-free evolution operator is Uo(t a ,tb) 
Texp ^— i / t * a dtH (t)^. These definitions ensure that at con- 
catenation level I the symmetrization procedure is applied to 
DD-free evolution segments of duration t; = T/G l , using a 
total of G l pulses. For piecewise constant Hamiltonians it 
has been shown that CDD achieves a = kh order decoupling 
P4ll48ll49l . Experimental studies support these conclusions 
El ED. 

QDDforAQC. — Another way to achieve high order decou- 
pling in a multi-qubit setting is QDD BTll . which is obtained 
by nesting two Uhrig DD [52 1 (UDD) sequences 



U^' M '\T,0) = U^ Q2 o U^ ni o Uo (T,0), (8a) 



t(M 2 ) 



(Mi) 



M+l 



fe=l 



where S { k M) = Tsm 2 [kn/(2M + 2)] and Q x * tl 2 e {X,Z} 
are the generators of S. The number of pulses in each se- 
quence, or sequence order, is {Mx,M2}, and dictates the de- 
coupling order a, requiring (Mi + 1)(M2 + 1) total pulses for 



a > min(Mi, M2) II53TI551 . In principle QDD uses exponen- 
tially fewer pulses than CDD to attain the same decoupling 
order, but as we shall see in fact CDD performs better in the 
context of our simulations of DD-protected AQC. 

Results. — For each algorithm (Grover, 2SAT) we used n = 
4 physical qubits to encode two logical qubits in the code 
C, and studied both CDD- and QDD-protected AQC. In our 
simulations the pulse intervals decrease as the total num- 
ber of pulses increases with the concatenation or QDD se- 
quence order, for each given value of the total time T. Start- 
ing for both algorithms from the uniform superposition state 
as the initial encoded state \i/j(0)), we computed \ip(T)) = 
U xDD (T,0) \ip(Q)}, where x=C or Q [Eqs. <§ and ([8a])], for 
a given Gaussian noise realization. Computing the ensemble 
average p(T) = (\ip(T))(ijj(T)\), we assessed performance 
using the trace-norm distance D[p(T), \<fr (T)}] J56), to 
quantify the difference between the encoded, xDD-protected 
state p(T) and the desired encoded final state |<i>o(T)), i.e., 
the ground state of H a d(T). In the ideal case of noise- and 
DD-free evolution |V>ad(*)) = Texp[-/ * H^t^dt'], the 
adiabatic theorem El guarantees D^tp^T)) , \$ (T))] « 
1 provided T »_ max s€{0) i } || ^H^/A^, s = t/T, 
In our simulations a finite range 

EE 



b e {1,2,3} 

of T's was used, and as can be seen from Figs 
D[p(T),\$ (T))] does indeed tend to zero for the ideal 
case as T is increased, though not monotonically in the 
Grover case llrjTI . The main effect of H m (t) is to cause 
D[p(T), |$o(T))] to diverge away from zero as T is in- 
creased, so that there is an optimal evolution time; this ef- 
fect of open system AQC, due to a competition between the 
benefit of increasing T for adiabaticity and the simultane- 
ous increasing damage due to system-bath coupling, has been 
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previously pointed out experimentally [9| and theoretically 
ifTZl l62l . The main role of DD protection, then, is to keep 
the fidelity of the AQC process as close as possible to the ideal 
and, in particular, to prevent min T6 [ ,T max ] D[p(T), l^o^))] 
from growing larger than some tolerance e > away from 
min T€[0jTmax] D[|^ ad (T)),|$o(^)>]- This is the sense in 
which DD-protected AQC approaches the ideal of FTAQC de- 
scribed in the introduction. 

CDD results for Graver's problem are shown in Fig. [TJ 
for increasing concatenation levels, at j3 = A/5. The faulty 
Grover evolution [generated by H^ d (t) + H err (t)] reaches a 
minimum deviation D « 0.13 at T « 4.5/A^ in and then di- 
verges from the ideal evolution [generated by H^(t)]. In 
contrast, CDD-protected evolution becomes remarkably close 
to the ideal evolution as the level of concatenation increases. 
Nearly ideal evolution is maintained essentially over the entire 
range of T values we simulated. 

Figure [2] shows our CDD results for the 2-SAT problem. 
The results are even better than those for the Grover prob- 
lem. CDD-protected evolution is essentially indistinguishable 
from the ideal at high enough concatenation level. The im- 
proved performance in the 2-SAT case can be attributed to the 
larger gap: A m ; n scales inversely with the number of qubits in 
2-SAT, instead of with the Hilbert space dimension as in the 
Grover case, rendering the 2-SAT evolution less sensitive to 
bath-induced noise effects. This is consistent with earlier ob- 
servations that algorithms associated with second order quan- 
tum phase transitions are more amenable to AQC than those 
for first order transitions lfT51 [32l [33l . 

Along with CDD, we analyzed QDD-protected AQC for 
Mi = M 2 = M e {1, 3, 6, 7, 14, 15}, where the odd sequence 
orders correspond to the first four levels of concatenation in 
CDD;, I e {1,2,3,4}, respectively, having the same number 
of pulses [(M + l) 2 and 4']. CDD and QDD are compared 
in our simulations at equal T values. QDDm performance 
improves monotonically for sufficiently large j3 (not shown), 
but we find that QDD15 performance is consistently inferior 
to that of CDD4 in the adiabatic regime of large TA, as can 
be seen in Fig. [3] This is surprising in light of the aforemen- 
tioned fact that in their non-AQC roles as quantum memory 
DD sequences, the respective decoupling orders of CDD; and 
QDDm are / and > M. We have confirmed that the results 
for the 2S AT problem are qualitatively similar, again favoring 
CDD. 

In Figs. [TJ and [2] the cutoff frequency j3 was fixed. Perfor- 
mance dependence on /3 is shown in Fig. [3] and is seen to be 
mild for CDD4, at each fixed value of T. The dependence 
on is further elucidated in Fig. |4j which shows that perfor- 
mance generally improves as f3 shrinks, as expected. In con- 
trast, QDD15 results show almost no dependence on /3, and 
the min r D[p(T(r)), |<&o (T(t) ) ) ] values are approximately 
twice as large as those obtained for CDD4 (not shown). In 
essence, the superior performance of CDD can be explained 
by recognizing that QDD is designed to suppress the system- 
bath coupling at the end of the pulse sequence, while CDD 
performs this suppression recursively all throughout the evo- 
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FIG. 3. (color online) CDD4 vs QDD15 for the Grover problem as 
a function of the normalized bath correlation time l/(/3T) and nor- 
malized total run time T/A° in . Increasing the bath correlation time 
1//3 at fixed T generally results in improved performance for CDD4. 
whose performance is significantly better than QDD15 in in the large 
T regime. All results were averaged over 30 realizations of H en (t). 
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FIG. 4. (color online) Performance of CDD4-protected evolution 
for the 2SAT problem, as a function of pulse interval r, for differ- 
ent values of the frequency cutoff f3 (r in units of 1//3 to separate 
the curves). Total time T = 4 4 r. The minimum in each fixed j3 
curve corresponds to an optimal pulse interval r opl (/3), and corre- 
sponding Top,(/3). Peak performance D[p(T opt (0)), |$o(T op ,(/3)))] 
improves as f3 is decreased, except at (3 = A^ in T /0.005 where self- 
averaging effects result from rapid fluctuations in t^(t) ("motional 
narrowing"). Results are averaged over 30 realizations of H al (t). 



lution. This is a better fit for AQC, with its time-dependent 
system Hamiltonian. 

Conclusion and future work. — We have introduced a high- 
order DD-based strategy for protected open-system AQC, and 
demonstrated using numerical simulations that it is capable 
of achieving high fidelities for a model of classical stochastic 
noise. At high enough concatenation level, the CDD-based 
protection strategy achieves fidelities which are essentially in- 
distinguishable from closed-system adiabatic evolution, for 
two algorithms associated with first and second order quan- 
tum phase transitions. Future work should elucidate the ques- 
tion of scaling of CDD resources with problem size, consider 
finite-width DD pulses, and generalize the results to a quan- 
tum bath. 
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